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Abstract 

Alzheimer's disease (AD) is a complex neurodegenerative disorder involving multiple cellular and molecular processes. The 
discovery of drug molecules capable of targeting multiple factors involved in AD pathogenesis would greatly facilitate in 
improving therapeutic strategies. The repositioning of existing non-toxic drugs could dramatically reduce the time and 
costs involved in developmental and clinical trial stages. In this study, preliminary screening of 140 FDA approved nervous 
system drugs by docking suggested the viability of the tricyclic group of antidepressants against three major AD targets, viz. 
Acetylcholinesterase (AChE), p-secretase (BACE-1), and amyloid (5 (AP) aggregation, with one member, protriptyline, 
showing highest inhibitory activity. Detailed biophysical assays, together with isothermal calorimetry, fluorescence 
quenching experiments, kinetic studies and atomic force microscopy established the strong inhibitory activity of 
protriptyline against all three major targets. The molecular basis of inhibition was supported with comprehensive molecular 
dynamics simulations. Further, the drug inhibited glycation induced amyloid aggregation, another important causal factor 
in AD progression. This study has led to the discovery of protriptyline as a potent multi target directed ligand and 
established its viability as a promising candidate for AD treatment. 
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Introduction 

Alzheimer's disease (AD) is the foremost cause of dementia in 
the ageing population affecting over 35 mUhon people worldwide. 
According to World Alzheimer Report 2013, this number is 
expected to increase by two fold in 2030. AD is a progressive 
neurodegenerative disorder that leads to the irreversible loss of 
neurons, intellectual abilities and eventually to death within a 
decade of diagnosis. Although the molecular bases of AD 
pathogenesis remains incompletely elucidated, the disease has 
been recognized as a multifactorial syndrome involving various 
molecular and cellular processes such as protein aggregation, 
oxidative stress, cell cycle deregulation and neuroinflammation 
[!]■ 

There are currendy several plausible hypotheses for AD 
pathogenesis. Cholinergic hypothesis states that the reduced 
cholinergic neurotransmission leads to the degeneration of 
cholinergic neurons and hence synaptic failure and cognitive 
dysfunction [2]. Following this, Acetylcholinesterase (AChE) was 
vahdated as a therapeutic target to reduce the degradation of 
acetylcholine in the synapse. AChE inhibitors (AChEIs) are 
effective in temporarily restoring cholinergic function, and 
constitute the majority of AD drugs currendy available in the 



market [3] . However, they are incapable of delaying or preventing 
neurodegeneration [4,5]. 

A variety of biochemical, genetic and pathological studies 
describe pivotal roles of the Amyloid P (AP) peptide in the 
pathogenesis of AD. "Amyloid hypothesis," describes the altered 
synthesis, aggregation and accumulation of AP which results in 
amyloid plaque formation [6]. While extracellular deposits of 
amyloid plaques are highly neurotoxic, recent studies have also 
implicated soluble, oligomeric aggregates of AP in neurotoxicity 
[7] . Therefore, a key therapeutic strategy for the treatment of AD 
involves the development of drugs targeted at inhibiting AP 
production, aggregation, destabilization and clearance of pre- 
formed fibrils [8] . A crucial step in AP production is the specific N- 
terminal enzymatic cleavage of the membrane embedded Amyloid 
Precursor Protein (APP) by the transmembrane aspartyl protease, 
P-secretase (BACE-1) [9-12]. Therefore, inhibiting BACE-1 has 
been considered as another attractive approach to prevent AP 
neurotoxicity. It is noteworthy that no inhibitors of AP aggregation 
or BACE-1 activity have reached the market yet, despite strong 
evidence of the causative roles of AP in AD. 

We further point out that AP is transported through a Receptor 
for Advanced Glycation End products (RAGE) and cause 
neuronal damage. Long-lived proteins are preferentially modified 
to form Advanced Glycation End products (AGE), and the stability 
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Figure 1 . Virtual screening by docl<ing. A. Heat map analysis of binding constants of 1 40 FDA approved nervous system drugs screened against 
A|3, AChE and p-secretase by Autodocl< tool 4.2. In the gradient ruler, red colour indicated strong binding (AG<— 6 kcal/mol), whWe green colour 
indicate weak binding (AG>— 3 kcal/mol) and the five drugs showing higher affinity to all the above mentioned targets were zoomed. B. Chemical 
structures of the five drugs. All are tricyclic anti-depressant drugs. 
doi:1 0.1 371 /journal.pone.01 051 96.g001 



of AP makes it an ideal substrate for non-enzymatic glycation and 
formation of AGEs [13]. In a recent study, it has been shown that 
AP-AGE formation may intensify the neurotoxicity whereas 
inhibition of this process significantly rescued the early cognitive 
deficit in mice [14]. Therefore, glycated AP has been considered as 
a more suitable ligand for RAGE, as it aggravates neuronal 
deterioration [14]. Hence inhibiting glycation of AP may be a 
valuable therapeutic strategy for AD. However, there has been no 
concerted effort to explore the inhibition of AP glycation as a 
therapeutic strategy. 

Efforts to target AChE inhibition [15]; AP production [16]; AP 
aggregation [17]; tau phosphorylation and aggregation [18] have 
been investigated largely in isolation, despite the complex nature 
of AD etiology. Recently, drug discovery in AD has gradually 
inclined towards development of "multi-target-directed ligands" 
(MTDLs) [19-22] which are efficient in treating complex diseases 
because of their ability to target multiple modes of disease 
pathogenesis. Further, to evaluate MTDLs for the AD treatment. 



"drug repositioning" seems to be an appealing strategy, as this 
approach has several advantages, including reduced time and cost 
necessary for clinical trials^'^ Priority candidate drugs for 
hypertension, retinoid therapy, diabetes and antibiotics with 
sufficient supporting evidences have been considered for reposi- 
tioning in AD [23]. However, to the best of our knowledge, 
repositioning drugs for multiple targets in AD is scarce. In this 
study, an in silico screening of 140 FDA approved drugs for 
neurological treatment was done against the primary targets of AD 
therapeutics, namely, acetylcholinesterase, BAGE-1, and AP 
aggregation. Further in vitro studies showed that amongst selected 
molecules, the tricyclic antidepressant protriptyline exhibited 
significant inhibition of AChE, as well as inhibition of other 
targets of AD. In addition, protriptyline was also found to inhibit 
glycation mediated AP aggregation. Mechanistic insights into 
protriptyline binding and inhibition of AChE, AP, and BACE-1 
activity were described in detail with molectolar dynamics 
simulation studies. 
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Figure 2. Protriptyline inliibits ACKiE activity. A. Determination of IC50 values of five drugs for AChE by using 0.05-0.8 miVl concentration range 
of all the drugs B. Estimation of the kinetic constants by Lineweaver-Burk analysis. AChE inhibition by protriptyline showed competitive inhibition. C. 
Isothermal Titration Calorimetric analysis of protriptyline - AChE interactions. The upper panel shows the raw data in the form of heat effect during 
titration and the lower panel shows corresponding thermogram representing the best fit curve D. Snapshot of drug binding with catalytic subsite of 
AChE E. snapshot of drug binding with anionic subsite of AChE F. Distribution of Protriptyline -anionic subsite (solid line) and Protriptyline -esteratic 
subsite (broken line) nonbonded (nonb) interaction energy; data are averaged over last 20 ns G. Fluorescence quenching of AChE by protriptyline H. 
CD spectra of binding of protriptyline to AChE and I. CD pro analysis to study the conformational change J. Evolution of the backbone RMSD for the 
Protriptylline bound (solid line) and free (broken line) AChE active sites from MD trajectories K. SASA distributions of active sites for Pro-bound (solid 
line) and free (broken line) AChE active sites from MD trajectories L. Measurement of AChE activity after treatment of neuro2a cells with 25 )j,M and 
60 |j.M protriptyline for 15 h. 
doi:10.1371/journal.pone.0105196.g002 
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Figure 3. Inhibition of Ap aggregation by protriptyline. APi3_22 aggregation in the absence and presence of protriptyline was investigated by 
A. Thioflavin T assay B. Light scattering C. CD analysis D. Atomic Force IVIicroscopy images [i and iv] of aggregated amyloid and Amyloid + 
protriptyline (0.010 mlVI) in 10x10 |jm^ and 20x20 |im^ surface area respectively. It is also represented in Line profile [ii and v] and 3D images [iii and 
vi]. 

doi:1 0.1 371 /journal.pone.01 051 96.g003 



Results and Discussion 

Tricyclic Antidepressant Drugs Display Strong Binding 
against Various Targets of AD In silico 

Multi-target-directed ligands (MTDLs) are likely to offer 
promising approaches for treatment of a disease as complex as 
AD [19-21]. The structures of 140 ligands were docked with the 
major targets of AD viz. AChE, BACE-1 and Afi aggregation. 
Ligands were scored based on electrostatic and hydrophobic 
contributions to the binding energy [24]. Furthermore, polar 



interactions were considered by H-bonding interactions analysis. 
Docking scores were used to rank ligands, depending on presence 
of number of H-bond donors and acceptors within the active sites 
[25]. Binding energy scores represented in the Heatmap 
(Figure lA) displayed variability in interactions of the ligands to 
the three targets of AD. There were several ligands that showed 
noteworthy interaction with at least two targets, but only few of 
them had strong interaction with all the targets. It was observed 
that five antidepressant drugs protriptyline, amytriptyline, mapro- 
tUine, doxepin and nortriptyline, which are tricyclic secondary 
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Figure 4. Destabilization of amyloid dimer by protriptyline. A. Evolution of monomer-monomer interaction strength over time for free dimer 
(brol<en line) and Protriptylline-bound dimer (solid line). Inset. Distributions of the interactions from multiple trajectories, and the dimer interactions 
with Protriptylline (in brown) B. Distributions of the asphericity for free (in broken line) and Protriptylline-bound (solid line) dimer C. Representative 
snapshot of most populated cluster of free, and D. Protriptylline-bound dimer [1 6-20 region in blue colour with 1 9-20 showed in line representation; 
protriptyline in red colour and two Ap peptides are in cyan and limon colour respectively] E. Residue-residue contact probabilities for free dimer, and 
F. Protriptylline-bound dimer G. Residue-wise Beta sheet percentages for free dimer (in red) and Protriptylline-bound dimer (in blue) H. Residue-wise 
helical percentages for free dimer (in red) and Pro-bound dimer (in blue). 
doi:1 0.1 371 /journal.pone.01 051 96.g004 
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Figure 5. BACE-1 inhibition by protriptyline. A. Determination of IC50 of BACE-1 by using various concentrations of protriptyline. The sigmoidal 
curve indicates the best fit for the percentage inhibition data obtained B. Lineweaver-Burl< analysis to estimate the kinetic constants. It showed 
competitive inhibition. C. Snapshot of drug binding with active site of BACE-1 . Active site residues in BACE-1 are in line representation. D. Active site 
of BACE-1. The structures from unbound (green) and ligand bound (orange) simulations are shown after all - atom superimposition. Snapshots are 
generated using PyMol. 
doi:1 0.1 371 /journal.pone.01 051 96.g005 



amines showed strong binding afFinity and broad specificity 
toward multiple targets of AD (Figure IB). 

Protriptyline Inhibits AChE Activity by Inducing 
Conformational Change in the Active Site 

AChE currently remains the foremost therapeutic target for AD 
and the current treatment and management of AD mainly involves 
use of acetylcholinesterase inhibitors [4,5] . Therefore, protripty- 
line, amitryptyline, maprotiline, doxepin and nortriptyline were 
assessed initially against AChE inhibition. While all the five drugs 
displayed inhibitory activity against this target (Figure 2A), 
protriptyline displayed highest inhibition, with the least IC50 of 
about 0.06 mM. In comparison, the other ligands displayed 
relatively higher ICsq values; the values corresponding to 
maprotUine, doxepin, nortriptyline and amitriptyline were 
0.1 mM, 0.480 mM, 0.135 mM and 0.6 mM, respectively. 
Hence, the inhibitory activity of protriptyline against the other 
targets of AD was evaluated in detail in the remaining study. 

The assessment of strong inhibitory activity of protriptyline was 
consolidated by enzyme kinetic studies, which suggested compet- 
itive inhibition of AChE by protriptyline. (Figure 2B). The 
apparent K,-,-, for AChE was determined to be ~0.025 mM by 
Lineweaver-Burk plot and it was found to be increased in the 
presence of protriptyline. The inhibition constant K; was 
determined from Cheng-PrusoflPs equation and found to be 



~0.001 mM. Thermodynamic studies of protriptyline-AChE 
interaction was carried out using Isothermal titration calorimetry 
(ITC) as it is one of the most widely used quantitative technique for 
direct measurement of the enthalpy change when two species 
interact, allowing the determination of heat of association, 
stoichiometry, and binding affmity from a single experiment 
[26,27]. The raw data and corresponding thermogram of the 
binding experiment is depicted in Figure 2C. Binding was 
strongly exothermic and showed 1:1 stoichiometry for AChE 
and protriptyline. The spontaneity of the process is evidenced by a 
negative change in the enthalpy, AH, and a positive change in the 
entropy, AS. Kinetic and ITC data demonstrated the binding of 
single molecule of protriptyline to the active site of AChE. The 
active site of AChE comprises two subsites; the anionic subsite 
(Trp82, Glul98, Tyr333) and the esteratic subsite (Serl99, Glu330 
and His443) [28]. Strong binding propensity of protriptyline to 
both the subsites of AChE was evidenced via the MD studies 
(Figures 2D and E). The mean interaction strengths of 
protriptyline with the anionic subsite is —31.3 (±9.1) kcal mol~ 
and —19.4 (±5.0) kcal mol with the esteratic subsite 
(Figure 2F). Breakup of the interactions into the non-bonded 
components showed that electrostatics plays a relatively stronger 
role than van der Waals interactions in protriptyline binding with 
active site. The mean interaction strengths of the active site 
residues with the ligand are provided in Table SI. 
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Figure 6. Protriptyline inKiibits glycation. A. Fluorescence emmission of AfJ and glycated Ap in presence of various concentration of 
protriptyline B. Light scattering and C. Kinetics of amyloid aggregation by Thioflavin T of A)3t3_22 and glycated Api3_22 in the absence and presence 
of protriptyline. 

doi:1 0.1 371 /journal.pone.01 051 96.g006 



Protriptylliiie-AChE interaction was also studied by fluores- 
cence spectrometry. A significant shift in tlie tryptophan fluores- 
cence emission spectra was observed upon protriptyhne binding 
suggested that the drug induces conformational change in AChE, 
thereby in its functional abilities. The fluorescence intensity of 
AChE exhibited an emission maximum at a wavelength, X^-^.^^, of 
325 nm in the unbound state. However, titration of the native 
enzyme with increasing concentrations of inhibitor resulted in a 
concentration-dependent quenching of the tryptophanyl fluores- 
cence (Figure 2G). A progressive red shift in the emission 
spectrum with increasing concentration of protriptyline suggested 
the conformation change in enzyme. Further, the conformational 
change was confirmed by CD analysis. CD analysis showed that 
the binding of protriptyline to AChE decreased the minima 
suggesting protriptyline induces conformational change in AChE. 
CDPro analysis revealed that binding of protriptyline increased P- 
sheets and decreased ot-helices in AChE. These results are very 
well in accordance with the fluorescence quenching experiment 
(Figure 2H and 21). 

The evidence of structural distortion was supported by the MD 
simulation data. In Figure 2J, we compare the backbone root 
mean squared deviation (RMSD) of the active site residues in the 
free and protriptyline bound state for a sample MD simulation 
trajectory. The higher RMSD in the latter is a distinct 
demonstration of the perturbative effect of the ligand on the 



active site structure. This structural distortion is supported by 
changes in the inter-residue distances of the active site residues (see 
Table S2). Distinct increases in several of the inter-residue 
distances are noted, especially in the distances involving Tyr333 
and His443. The sharpest increase in the mean inter-residue 
distance is found for two residues belong to the esteratic subsite, 
namely Ser 1 99 and His443 that play critical roles in the hydrolysis 
of acetylcholine. The increase from 7.8 A in the free state to 
1 1 .5 A in the bound state strongly suggests that protriptyline 
binding directly affects the catalytic ability of AChE. Further, we 
found protriptyline binding at the active site to be commensurate 
with an increase in its solvent accessibility. The solvent accessible 
surface area, or SASA, was calculated by running a spherical 
probe of 1 .8 A radius around the surface and calculating the area 
covered by the probe / calculated in the standard manner. In 
Figure 2K, we have presented histograms of the SASA obtained 
from the MD trajectories of the free and protriptyline bound 
AChE. The active site SASA increased from 1533 (±38.9) to 
1623 (±33.0) A^ upon protriptyline binding. 

Additionally, to confirm the in vitro results, the effect of 
protriptyline on AChE activity was assayed in neuro2a cell line. 
Protriptyhne treated cells showed concentration dependent 
decrease in AChE activity (Figure 2L). These results demon- 
strated that the in vilro IC50 of protriptyline for AChE was 
sufficient to inhibit its activity in cell culture. Thus, this study 
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doi:10.1371/journal.pone.0105196.g007 
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Figure 8. Cell viability in neuro2a cells. Effect of various concentrations of protriptyline (25-500 nlVI) on cell viability was assessed by MTT assay. 

Cells were 90% viable up to 150 jiM protriptyline concentration. 

doi:10.1371/journal.pone.0105196.g008 



suggested that protriptyline inhibits AChE activity in vitro as well 
as in neuro 2a cells, by binding to the active site and causing 
conformational change. 

Protriptyline Inhibits A/? Self-Assembly 

Recent studies suggest that the KigLVFFjo segment in the AP 
sequence is crucial for the peptide's oUgomeric properties as well 
as fibriUogenetic behavior [29-33]. The sequence 
HHQKLVFFAE corresponding to APi3_22 was used for our in 
vitro amyloid aggregation inhibition studies. P-sheet rich structures 
are a common feature of amyloid aggregates [34] and bind to the 
molecule Thioflavin T (ThT). Therefore, ThT fluorescence assays 
are frequently used to monitor aggregation of amyloidogenic 
peptides [35]. In Figure 3A, we present ThT fluorescence as a 
function of time for APi 3_22 in the absence or presence of 10 |J,M 
protriptyline. It is obs(T\cd that the lag time corresponding to the 
nucleation phase is increased from ~25 to about ~50 hours in the 
presence of protriptyline. The data further shows a dramatic 
inhibition in the growth phase in the presence of protriptyline. The 
fluorescence intensity corresponding to the saturation phase was 
further lowered in the presence of protriptyline, demonstrating the 
ability of protriptyline to reduce AP aggregation. Concentration 
dependent decrease in Thioflavin T fluorescence was observed in 
presence of protriptyline at 7* day (Figure SI). Static light 
scattering experiments were performed to investigate the relative 
decrease in average molecular mass of AP aggregates in the 
presence of protriptyhne (Figure 3B). As expected, protriptyline 
reduced the light scattering in a concentration dependent manner. 
Furthermore, Far-UV CD spectra of the APi3_22 were recorded 
with and without protriptyline to monitor any possible alterations 
in secondary structural propensities (Figure 3C). CDPro analysis 
comparing the spectra obtained at the two conditions showed that 
protriptyline treatment reduced P-sheet formation and increased 
overall helicity. The reduced P-sheet formation by CD analysis 
corroborated the Thioflavin T results. We further performed 
atomic force microscopy (AFM) experiments on the peptide 
aggregates formed with and without protriptyline treatment. The 
results of an AFM scans for APi3_22 incubated for 7 days with and 
without protriptyhne over scanning areas of 10 |J,mxlO |im and 
20 |J,mx20 |J,m are compared in Figure 3D. 

First, the AFM profiles suggested a distinct reduction in the 
fibrillar density and a high degree of size dispersion resulting from 
protriptyline treatment. Further, the line profile and 3D AFM 
images showed that the average height of the APi3_22 aggregates 
was ~400 nm, with the maximum height reaching ~600 nm. The 



height of the aggregates reduced to ~ 1 00 nm in case of 
protriptyline treated samples. It was quite convincing from the 
AFM images that the amyloid protein is aggregated and 
protriptyline treated sample protein showed decreased aggrega- 
tion. 

For a m()let:ular level insight into the inhibitory action of 
protriptyline on AP self-assembly, we performed MD simulations 
of AP dimerisation in the absence and presence of protriptyline 
molecules. A P-sheet rich full-length monomer was considered for 
these studies. A single protriptyline molecule was found to bind to 
the AP monomer with a mean binding strength of —84.0 (±34.0) 
kc:al moF' (Figure S2A). The mean monomer-monomer 
interaction strength obtained at the end of multiple independent 
trajectories over a combined total simulation time of 60 ns was — 
249.1 (±87.1) kcal mol~'. In the presence of protriptyline, 
however, the monomer-monomer interaction weakened signifi- 
cantiy, with a mean value of only —69.7 (±31.5) kcal mol In 
Figure 4A, we have depicted evolution of monomer-monomer 
interaction over simulation time for sample simulations in the 
absence and in the presence of protriptyline; the distributions of 
these interactions, along with distribution of the protriptyline- 
dimer interaction is depicted as an inset. The weakening of the 
inter-monomer interactions due to protriptyline binding results in 
the formation of complexes that are less compact compared to the 
pure dimeric forms. In Figure 4B, we compare distributions of 
asphericities ly, of the pure and protriptyline-bound dimer 
complexes; ly value of 0.0 denotes perfect sphericity while 
increasing values denote increased asphericity. The protriptyline 
bound dimer complexes have a wider distribution and a higher 
mean value of ly; the mean asphericity values of the free and 
protriptyline -bound complexes are 0.2 and 0.5, respectively. 
Further, the mean radii of gyration (Rg) of these complexes were 
found to be 13.0 (±0.1) and 16.5 (1.0) A, respectively. The 
difference in dimeric compactness upon protriptyline binding is 
evident from representative snapshots depicted in Figures 4C 
and 4D. Additionally, clustering analysis was done to identify the 
key residues of AP interacting with the protriptyline during 
molecular dynamics simulation. A representative snapshot from 
most populated cluster is shown in Figure 4D. In two of the most 
populated clusters (70% of total snapshots used in the clustering 
analysis), protriptyline interacts most strongly with the KLVFF 
region. 

In order to evaluate the effects of protriptyline binding on the 
inter-monomer associations, we have analyzed the nature of inter- 
residue contacts. The average numbers of inter-residue side-chain 
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Figure 9. Protriptyline as MTDL. The scheme represents that protriptyline (at the center) Is able to Inhibit key targets of AD pathogenesis such as 
AChE, BACE-1, Amyloid aggregation and glycatlon Induced amyloid aggregation. 
dol:1 0.1 371 /journal.pone.01 051 96.g009 



cross-contacts were reduced up to 46% upon protriptyline 
binding. Here, as in previous studies [36], we have defined two 
residues to be in contact if the maximum separating distance of 
their side-chains does not exceed 7 A. In Figures 4E and 4F, we 
compare the inter-residue contact probabilities for tlie free and 
protriptyline-bound complexes from the MD data. The largest 
numbers of contacts between the KLVFF regions are lost upon 
protriptyline binding. However, we note the emergence of a small 
extent of extraneous contacts, particularly involving the N- and C- 
terminal regions of the different monomeric units. 

Further the effect on the secondary structural propensities of the 
A|3 units due to protriptyline binding was investigated. Protripty- 
line was found to reduce (3-sheet and induce higher helical 
propensities in the monomeric form of AP (see Figure S2 B). 
Figure 4G and 4H present residue-wise P-sheet and helical 
propensities of the free and protriptyhne-bound dimeric complex- 
es, respectively, from the MD data. A sharp decrease in /?-sheet 
propensity is found uniformly along the AP sequence, including in 
the residue span H13HQKLVFFAE22. The decrease in P-sheet 
propensity is accompanied with an overall increase in helical 
conformations. Sharp increase in helicity was observed near the N- 
terminal and KLVFF regions. The alterations to secondary 
structural propensities thus observed from MD analysis are an 
excellent corroboration of the CD and ThT binding data. 



Protriptyline Inhibits BACE-1 Activity 

BACE-1 is a key enzyme required for AP production. Hence, 
BACE-1 inhibition is an attractive target for countering AD [10- 
12]. BACE-1 assay demonstrated decreased activity with increas- 
ing concentrations of protriptyline having IC50 ~0.025 mM 
(Figure 5A). Protriptyline inhibited BACE-1 by competitive 
inhibition as depicted by Lineweaver-Burk plot (Figure 5B), 
Apparent K,,-, and K; of BACE-1 were calculated as mentioned 
above for AChE and found to be 0.0025 mM and 0.005 mM 
respectively. Apparent K.^ of BACE-1 was increased in the 
presence of protriptyline Competitive inhibition was evidenced by 
MD simulation analysis that illustrated protriptyline binds strongly 
at the active site of BACE-1 comprised of Asp32 and Asp228 
(Figure 5C). The mean binding strength required for protripty- 
line to bind to the active site is —29.5 (±7.0) kcal mol '. The 
center of mass distance between the two residues increases from 
5 A to 8 A as a result of protriptyline binding (Figure 5D). The 
comparison of the root mean squared deviation (RMSD) of C^ 
atoms of the active site in the unbound state with the protriptyline 
bound state indicated structural distortion of the arrangement of 
the active site (see Figure S3A). The binding was also found to 
induce significant alterations to the local secondary structural 
propensity around the active site (see Figure S3B). Therefore, 
inhibition of BACE-1 by protriptyline is an extra benefit as it 
prevents AP generation. And even if there is some production, 
protriptyline will obstruct it to get aggregate. 
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Protriptyline Inhibits Glycation Associated Aggregation 
of AP 

AD is also referred as type III diabetes [37] and its pathogenesis 
has been correlated with the extent of glycation [38]. Recent 
studies suggested glycated AP is more neurotoxic than native AP 
[13]; therefore the effect of protriptyline on glycation of AP was 
investigated. Glycated proteins emit fluorescence at 440 nm upon 
excitation at 370 nm. Fluorescence assay illustrated that AP 
undergoes glycation. The increase in glycation associated fluores- 
cence was reduced by protriptyline in a concentration dependent 
manner (Figure 6A). Glycation enhances the aggregation and 
also alters the secondary structure of proteins. Static light 
scattering was used to study glycation induced protein aggregation. 
Concentration dependent decrease in light scattering was observed 
(Figure 6B). Further, it was studied by Thioflavin T fluorescence 
assay. Kinetics of aggregation displayed increased Thioflavin T 
fluorescence during glycation reaction. Lag phase for AP and 
glycated AP aggregation was increased in the presence of 
protriptyline with decreased Thioflavin T fluorescence (Fig- 
ure 6C). In addition to inhibition of AP glycation, the drug also 
inhibited glycation of insulin and BSA(Bovine Serum Albumin) 
evidenced by MALDI, AGE fluorescence and Thioflavin T 
fluorescence assay (See Figure S4 and S5). 

Protriptyline Does Not Affect Other Proteases 

Protriptyline inhibited multiple targets of AD; therefore we 
further studied the influence of this drug on other enzymes such as 
trypsin and Qt-secretase. It was interesting to observe that it was not 
able to inhibit trypsin activity even at 0.1 and 0.5 mM 
(Figure 7A). For ADAM 17, protriptyline was found to be a 
weak inhibitor. It was showing ~3% inhibition at 0.1 mM and 
~25'yo inhibition at 0.5 mM (Figure 7B). It suggested that 
protriptyline is not a non-specific inhibitor as it was not able to 
inhibit other proteases. 

Viability of neuro2a cells in the presence of protriptyline 

The effect of various concentrations of protriptyline (25- 
500 nM) on cell viability was evaluated by MTT assay in neuro2a 
neuroblastoma cells. It was observed that there was more than 
90% of cells were viable upto 150 jiM of protriptyline, and a 
drastic reduction in cell viability was observed at 200-500 ^lM 
protriptyline concentrations (Figure 8). In this study the IC50 of 
protriptyline for all the three targets AChE, AP, BACEl was less 
than 150 |iM, this concentration was found to be non toxic to 
cells. 

Conclusion 

The multifactorial nature of AD makes its treatment complex 
and unmanageable. The discovery of molecules that can inhibit 
multiple pathways of the disease should significantly advance 
therapeutic strategies. In this study, we investigated the efficacy of 
the tricyclic antidepressant, protriptyhne, against important AD 
targets. Our in vitro and in silico investigations established the 
inhibitory effects of the drug on AChE, amyloid aggregation, 
BACE-1 and glycation (Figure 9). Protriptyline was able to 
inhibit AChE and P-secretase by binding at the active site and 
causing conformational changes. In addition, it strongly prevented 
self-assembly of AP and glycated Ap. It is a FDA approved drug 
for the treatment of depression, narcolepsy. Attention Deficit 
Hyperactivity Disorder (ADHD) and headaches and its ability to 
cross blood brain barrier (BBB) [39] is an additional advantage, 
which is a crucial requirement of molecules used for intra-cranial 
diseases. Furthermore, as there is high prevalence rate (30-50%) of 



AD and depression co-morbidity, the use of antidepressants can be 
a rational complementary therapy for AD treatment [40]. 
Therefore, antidepressant activity of this drug could be an added 
advantage when dealing with AD complications related to 
depression. To the best of our knowledge, this is the first study 
in which an anti-depressant drug has been shown to inhibit 
multiple targets of AD. Our results strongly ratify protriptyUne as a 
promising candidate for AD therapy, and its further evaluation in 
animal and cUnical studies. 

Experimental Section 

Materials 

AU chemicals were procured from sigma unless and otherwise 
stated. 

Virtual Screening for Multi-target Drug Ligand for AD 

Structures of 140 FDA approved antiepileptics, psycoleptics, 
analgesics, psycoanaleptics, anti-Parkinson and other nervous 
system drugs were obtained from DrugBank (The DrugBank 
database website. Available: http://www.drugbank.ca/. Accessed 
2013 Feb 19) database and optimized for their 3D coordinates 
using Marvin Sketch Tool (CheniAxon website. Available: http:/ / 
www.chemaxon.com. Accessed 2013 Feb 21). Three dimensional 
structures of Human acetylcholinesterse (PDB ID: 1B41), P- 
secretase (PDB ID: 2HM1) and AP peptide (PDB ID: IZOQ) were 
accessed from RCSB PDB. Protein structures were energy 
minimized using Swiss PDB viewer (Swiss-Pdb Viewer application. 
Available: http://spdbv.vital-it.ch/. Accessed 2013 Feb 21). 
AutoDock tool 4.2 [41] was used to conv(;rt nx eptor and ligand 
from *.pdb to *.pdbqt format and to set other docking parameters. 
Grid map was set around the active site of acetylchoUnesterase 
(Serl99, Glu330 andHis443), P-secretase (Asp32 and Asp228) and 
KLVFF (residue 17 to 21) region of AP protein involved in 
aggregation. Virtual screening was carried out using AutoDock 
Vina software and the Lamarckian genetic algorithm as a 
searching procedure [42]. Binding energy obtained for each 
complex was represented in heat map format using MeV software 
packages (MeV: MultiExperiment Viewer. Available: http:/ /www. 
tm4.org/mev/. Accessed 2013 Mar 5) [43]. The gradient ruler is 
an indicator of interaction strength. Molecules showing strong 
binding against all the selected targets were selected for further in 
vitro and molecular simulation studies. 

Acetylcholinesterase Inhibition Assay 

The modified method of EUman et al. [44] was adopted to 
measure AChE activity. Briefly, 25 |J,1 (0.3 U/ml) AChE from 
Electric eel fish {Electrophorus electricus) was incubated with and 
without different concentrations of drugs selected from Molecular 
Docking studies. The reaction was carried out for 15 min at 25°C. 
500 |xl of 5, 5-dithiobis (2-nitrobenzioc) acid DTNB (3 mM) was 
then added and reaction was initiated by the addition of 100 Jtl 
substrate acetyl thiocholine iodide (ATCI) (15 mM). Total volume 
of the reaction was made up to 1 ml by Tris buffer, pH 8.0. ATCI 
hydrolysis was measured by colored product formation, 5-thio-2- 
nitrobenzoate anion by reaction between DTNB and thiocholine, 
a hydrolysis product of ATCI. The formation of the colored 
product was measured at 405 nm wavelength after 30 min. The 
background was determined from negative controls (omission of 
AChE enzyme). 

Isothermal Titration Calorimetry 

ITC was performed using Microcal Auto-iTC instrument (GE 
Healthcare). 40 injections of 2 (il protriptyUne (Stock = 2.2 mM) 
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was titrated against 0.3 U/|j.l solution of AQiE. Experiments were 
carried out at 25°C in a Tris buffer, pH 8.0. Reference titration 
was carried out by injecting the same concentration of protrip- 
tyline into buffer. Reference titration was subtracted from 
experimental titration. Origin 6.0 software was used to derive 
affinity constants (Kj), the molar reaction enthalpy (AH) and the 
stoichiometry of binding (N), by fitting the integrated titration 
peaks. 

Fluorescence Analysis of AchE-Protriptyline Interaction 

AchE-protriptyline interaction was also studi(;(l b}' measuring 
tryptophan fluorescence using Varioscan plate reader. AChE was 
excited at 280 nm and emission was scanned from 300 nm to 
500 nm. Titration of enzyme with protriptyline was performed by 
the addition of different concentrations of inhibitor (10 )iM- 
100 nM) to a fixed concentration (0.2 U/|j.l) of enzyme solution. 
Background buffer and inhibitor spectra were subtracted and 
graphs were smoothed. 

BACE-1 Inhibition Assay 

BACE-1 activity was studied in accordance with the manufac- 
turer's instructions (Sigma). To test the effect of protriptyline on 
BACE-1 activity, 1.8 U of BACE-1 enzyme and 10 |a,M of |3- 

secretase specific peptides conjugated to fluorogenic reporter 
molecules were incubated with or without various [:oncentrations 
of protriptyline for 2 h at 37"C. Negative control included all the 
reactants except BACE-1 enzyme. After 2 h, fluorescence emission 
was measured at 405 nm upon excitation at 320 nm. 

Inhibition Kinetics for AChE and BACE-1 

Michaelis-Menten constant (K„) was determined by measuring 
the activities of AChE and BACE-1 using various concentrations 
of ATCI (100 nM-1000 nM) and BACE-1 substrate (0.25 nM~ 
6 H-M) respectively. Lineweaver-Burk double reciprocal plot was 
plotted in order to determine the K^. The protriptyline inhibition 
kinetics was analyzed over a range of concentration (50 |J.M- 
1500 |XM) and (10 |J,M-150 |xM) for AChE and BACE-1 respec- 
tively. The IC50 of protriptyline for both the enzymes was 
calculated by determining the inhibitor concentration at which the 
enzyme activity is 50% inhibited. The K; was calculated directiy 
from IC50 value using Cheng-Prussoffs classical equation [45]. 

In order to determine the type of inhibition, AChE and BACE- 
1 were incubated with 60 |iM and 15 |J,M protriptyline concen- 
trations respectively and assayed at increasing concentrations of 
AChE substrate (ATCI, 100 |xM-1000 \iM) and BACE-1 sub- 
strate (1 |J,M-10 |J.M) respectively. The reciprocals of reaction rate 
(1/v) for each inhibitor concentration were plotted against the 
reciprocals of the substrate concentrations (1/S). Mode of 
inhibition by protriptyline was determined from the graphical 
representation. 

Inhibition of Ap Aggregation 

13-22 peptide (HHQKLVFFAE), the aggregation prone 
region of AP, was synthesized from Thermo Fisher Scientific. 
Synthetic APi3_22 peptide was dissolved in 10% ammonium 
hydroxide and sonicated for 5 min. It was then diluted in 10 mM 
PBS, pH 7.0 to a final concentration of 200 ^M. 100 |J,M of APi3_22 
was incubated with and without various concentrations of protrip- 



tyline in 10 mM PBS, pH 7.0 at 37 C for a week and these samples 
were further used for aggregation inhibition assays. 

Glycation Inhibition Assay 

Anti-glycation activity of protriptyline was studied using insulin, 
BSA and Ap. The details of insulin and BSA glycation are described 
in Method SI. In Ap glycation studies, 100 jlMAPi 3.22 and 0.1 M 
glucose were incubated with and without various concentrations of 

protriptyline in 10 mM PBS, pH 7.0 for 7 days. These samples were 
further used to study glycation mediated AP aggregation by 
different assays. 

Thioflavin T Assay 

Thioflavin T assay was performed on 7* day of incubation. 
50 |il (25 |xM) of APi3_22 was mixed with 150 |il of Thioflavin T 
(ThT) stock solution (50 ThT in PBS pH 7.0) and placed in 
96-weIl plate (black with flat bottom. Comings). Fluorescence 
emission was measured at 460-550 nm upon excitation at 
440 nm. To account for background fluorescence, the fluores- 
cence intensity measured from each control solution without AP 
was subtracted from each solution containing APi3_22- Similarly, 
Thoflavin T assay was performed for aggregation kinetics of APi 3_22 
and APi3_22 glycation with and without protriptyline. 

Light Scattering Analysis 

Inhibition of aggregation/ glycation mediated aggregation of 
Ap was detected by static light scattering method using a Perkin- 
Elmer Luminescence spectrometer LS50B. Both excitation and 
emission wavelengths were set at 400 nm. Excitation and emission 
slit width was set to 10 nm and 2.5 nm, respectively. Scattering 
was recorded for 60 sec. 

Circular Dichroism Spectroscopy 

The far UV CD spectra (in wavelength range of 190-250 nm) 
of AP 13-22 (20 Hg/ml) with and without protriptyline was 
recorded on a Jasco-J815 spectropolarimeter at ambient temper- 
ature. In case of AChE, the enzyme was incubated with 
protriptyline and CD spectra were recorded at 0 min and 
30 min. Each CD spectrum was accumulated from thr(^e scans 
at 50 nm/ min with cell path length of 0. 1 cm. Contribution due to 
buffer was corrected in all spectra and observed values were 
converted to mean residual ellipticity (MRE) in deg cm^ dmol~ 
defined as 

MRE= -rr—r- 

lOdcr 

Where M is the molecular weight of the protein, 9x is CD in 
mUlidegree, d is the path length in cm, c is the protein 
concentration in mg/ml and r is the number of amino acid 
residues in the protein. Secondary structure content of the amyloid 
with and without protriptyline was calculated using the CDPro 
software (CDPro software package. Available: http://lamar. 
colostate.edu/~sreeram/CDPro/main.html. Accessed 2013 Oct 
18). 

Atomic Force Microscopy 

For atomic force microscopy (AFM) analysis, 10 |il of each 
sample was deposited on a piece of freshly cleaved mica disk. The 
disk was washed with water and dried overnight. The sample was 
mounted onto a Multimode scanning probe microscope equipped 
with a Nanoscope IV controller from Veeco Instrument Inc., 
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Santa Barbara, CA. All tJie AFM measurements were done under 
ambient conditions using the tapping-mode AFM probes model - 
Tapl90Al purchased from Budget Sensors. The radii of tips used 
in this study were less than 10 nm, and their height was ~17 |j,m. 
The cantilever used had a resonant frequency of ca. 162 kHz and 
nominal spring constant of ca. 48 N/m with a 30 nm thick 
aluminium reflex coating on the back side of the cantilever of the 
length 225 [im. For each sample, three locations with a surface 
area of 20x20 l^m^ and 10x10 \lm^ for amyloid and protriptyline 
treated amyloid were imaged with a frequency of 1 Hz and at a 
resolution of ,5I2x,512 dpi. Representative images were selected 
for comparative studies. 

Measurement of Glycation Associated Fluorescence 

Glycation associated fluorescence of AP was measured in 
amyloid, glycated amyloid treated with or without protryptUine at 
370 nm excitation and emission was scanned from 400-550 nm. 

BApNA Assay 

Activity oi Bovine trypsin was estimated using enzyme-specific 
chromogc-nic substrate BA/*NA [46]. In brief, 10 |ig Bovine trypsin 
was incubated with and without 100 ^M, 200 jiM, 500 nM 
protriptyline at 37°C for 15 min and volume was made upto 
150 Ml with 0.1 M Tris-HCl pH 7.8. Further, 1 ml BA/)NA was 
added to the reaction mixture and incubated for 10 min at 37°C. 
The reaction was terminated by addition of 200 \il of 30% acetic 
acid and absorbance was measured at 410 nm. 

ADAM 17 Assay 

ADAM 1 7 assay kit (Enzo Life Sciences) was used to study the 
effect of protriptyline on ADAM 17. Briefly, ADAM 17 and 
fluorogenic peptide (substrate) was incubated with and without 
100 nM protriptyline at 37°C for 10 min. The negative controls 
had the reaction mixture without enzyme. Fluorescence was 
measured at 328 nm and 420 nm for excitation and emission 
respectively. 

Cell culture 

Murine neuro2a neuroblastoma cells were obtained from 
National Center for Cell Science, Pune, India. The cells were 
maintained in Dulbecco's modified Eagle's medium (DMEM) and 
10% Fetal Bovine Serum (FBS). Cells were maintained at 37°C in 
humidified air containing 5% CO2 and were grown in monolayer 
cultures. 

Determination of AChE inhibition in cultured cells 

Equal numbers of cells were seeded in 25 cm^ flasks (1.5 
million). Cells were allowed to adhere and attain their morphol- 
ogy. Cells were serum starved for 24 h and treated with or without 
25 |a,M or 60 |a,M of protriptyline concentration for 16 h. Cells 
were trypsinised and given two washes with ice cold PBS. Protein 
was extracted by sonicating the cells for 20 min. The cell lysate 
was ccntrifuged for 60 min at 16000 rpm. Supernatant was 
collected and protein concentration was determined by Bradford's 
method. Acetylcholinesterase activity was assayed as described 
above by Ellman's assay. 

Cell viability 

Cell viabiUty following exposure to protriptyline was measured 
by MTT reduction assay. Neuro 2a (N2a) neuroblastoma cells 
were seeded at a cell density of 1 x 10* cells per well in a 96 well 
plate. After the cells adhered and attained their morphology, they 
were serum starved for 24 h prior to treatment with various 



concentrations (25-500 [iM) of protriptyline in triplicate for 15 h. 
After incubation, cells were given one wash with PBS and 100 jll 
fresh serum free media was added. 20 |a,l of 5 mg/ml MTT 
(dissolved in PBS) was added to each well and incubated in dark at 
37^0 until violet formazan crystals were observed. Media from 
each well was discarded and crystals were dissolved in 100 Jtl 
DMSO. Absorbance was measured at 555 nm using Biorad iMark 
microplate reader. 

Molecular Dynamics Simulations 

All simulations in this study were carried out with the 
NAMD2.9 package [47], using the CHARMM22 all-atom force 
field with CMAP correction for the proteins [48,49]. Force field 
parameters for protriptyline were generated using the Swis- 
sPARAM tool [50], and refined via electronic structure calcula- 
tions using Gaussian03. This strategy has been used in several 
recent studies [51-53]. Simulations were carried using a time step 
of 2 fs in the isothermal-isobaric (NPT) ensemble at a temperature 
of 310 K and a pressure of 1 atmosphere. Each system was 
sampled for a total duration of 60 ns with multiple trajectories. 
The SHAKE algorithm [54] was used to constrain bond lengths 
involving hydrogen atoms. Constant temperature was maintained 
with Langevin dynamics with a coUision frequency of 1 ps and 
constant pressure was maintained using the Langevin piston Nose- 
Hoover method [55]. Three-dimensional orthorhombic periodic 
boundary conditions were employed and full electrostatics 
calculated with the particle-mesh Ewald method [56]. A non- 
bonded cutoff distance of 12 A was employed, which were 
smoothened at a distance of 10.5 A. Details of system setups and 
trajectory analysis are provided in Method S2. 

Statistical Analysis 

All the experiments performed independently three times. 
Student's t-test was used for statistical analysis. Data were 
expressed as mean ±SD. A p-vaiue <0.05 was considered as 
statistically significant. 

Supporting Information 

Figure SI Thioflavin T Assay. Protriptyline causes concen- 
tration dependent decrease in A(3 aggregation. 
(DOCX) 

Figure S2 A. Snapshot of drug binding with Ap mono- 
mer. Central hydrophobic core are in line representation B. 
Comparison of residue wise percentage of helix from 
unbound and ligand bound Ap monomer simulated 
trajectory. 
(DOCX) 

Figure S3 A. Root mean square deviation (RMSD) of 
active site regions of P-secretase in unbound (red) and 
ligand bound (blue) simulated trajectory. B. Compari- 
son of residue wrise percentage of beta from unbound 
(red) and ligand bound (blue) simulated trajectory. 
(DOCX) 

Figure S4 A. Control insulin B. Glycated insulin C. 
glycation inhibition in presence of 500 [lM. and D. 
1000 |iM protriptyline. These spectra were acquired on a 
positive reflector mode by MALDI-TOF-MS. Glycated peaks are 
shown by black arrow. 
(DOCX) 

Figure S3 BSA (Bovine Serum Albumin) glycation inhi- 
bition assay. Glycation inhibition of BSA was studied by A. 
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Measurement of AGE fluorescence, % glycation inhibition was 
plotted and B. Thioflavin T fluorescence assay, concentration 
dependent decrease in thioflavin T fluorescence was observed. 
(DOCX) 

Table SI Acetylcholinesterase-protriptyline interaction 
energy calculations. Interaction energy of ligand with active 
site residues averaged over last 20 ns of all simulated trajectories. 
Interaction strengths are in kcal mol-1 unit. Standard deviations 
are provided within braces. 
(DOCX) 

Table S2 Inter-residue distance between active site of 
AChE residues averaged over last 20 ns of free (vertical) 
and ligand bound (horizontal) simulated trajectories. 

Distances are in A unit. Standard deviations are provided within 

braces. 
(DOCX) 

Method SI Inhibition of Insulin and BSA glycation. 

(DOCX) 
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